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ABSTRACT 



In polarimetry it is important to characterize the polarization properties of the instrument itself to 
disentangle real astrophysical signals from instrumental effects. This article deals with the accurate mea- 
surement and modeling of the polarization properties of real aluminum mirrors, as used in astronomical 
telescopes. Main goals are the characterization of the aluminum oxide layer thickness at different times 
after evaporation and its influence on the polarization properties of the mirror The full polarization prop- 
erties of an aluminum mirror are measured with Mueller matrix ellipsometry at different incidence angles 
and wavelengths. The best fit of theoretical Mueller matrices to all measurements simultaneously is ob- 
tained by taking into account a model of bulk aluminum with a thin aluminum oxide film on top of it. Full 
Mueller matrix measurements of a mirror are obtained with an absolute accuracy of ~ 1 % after calibration. 
The determined layer thicknesses indicate logarithmic growth in the first few hours after evaporation, but 
it remains stable at a value of 4.12±0.08 nm on the long term. Although the aluminum oxide layer is 
established to be thin, it is necessary to consider it to accurately describe the mirror's polarization proper- 
ties. 

Subject headings: Astronomical Instrumentation 



1. Introduction 



Polarimetry is a powerful astronomical technique for characterizing e.g. magnetic fields and scattering media. 
The careful analysis of the polarization of light is, however, hampered by the instruments themselves as individual 
optical components induce and modify polarization. To improve the polarimetric accuracy of current and future po- 
larimeters, it is therefore crucial to know and understand the polarization properties of the telescope and the instrument. 
All reflecting surfaces such as the aluminized telescope mirrors produce some instrumental polarization at oblique in- 
cidence angles. The precise polarization properties of aluminum mirrors are notoriously difficult to model, because of 
variations of the refracti ve index, the presence of an aluminura oxid e (AI2O3) film on the mirror, and the presence of 
dust on the surface (e.g.l Harrington & Kuhn 2008 : Joos et al. 2008 ). This article deals with the empirical character- 
ization and physical modeling of the non-ideal behavior of aluminum mirrors due to the aluminum oxide layer. The 
second paper in this series discusses the influence of dust on the surface (Snik et al. 2009; in preparation). 



The importance of taking into account the aluminum oxi de layer was theoretically shown b 



e.g.. lBurge & Bennett 

(119641) . Ellipsometry with a He-Ne laser was performed bv lSankarasubramanian et al.l (119991) . who showed the im- 
provement of the Mueller matrix model after taking into account the aluminum oxide layer. However, the deduced 
aluminum oxide thickness values (~ 50 nm) are an order of magnitude larger than what is expected from a layer formed 
by qua ntum mechanical tunneling of electrons through the layer at room temperature. This theory, developed bv lMott 
(I1939I) . predicts logarithmic layer growth in the first few days, up to a thickness of about 5 nm. This theoretical model 



-2- 



successfuUy explained the thickness measurements performed bv lJeurgens et alj (l2002h based on electron microscopy 
and X-ray photoelectron spectroscopy, which yields an amorphous aluminum oxide layer with a uniform thickness 
of 0.5-4.0 nm for low temperatures (< 573 K). To fully characterize the aluminum oxide layer and its influence on 
an aluminum mirror's polarization properties as a function of wavelength and time after evaporation, we constructed 
a precision ellipsome ter based on a broad-band light s ource in combination with interference filters. Compared to 



Sankarasubramanian et al. (Il999 ^) at a single wavelength and a single time after evaporation. 



the measurements of 

we studied the growth of the oxide layer, by accurately measuring a mirror's Mueller matrices at different incidence 
angles, at several wavelengths, and at different points in time after evaporation. 



2. Measurements 

A complete Mueller matrix ellipsometer, based on liquid crystal variable retarders (LCVRs), was developed for 
use in transmission, reflection and scattering measurements. Any incidence angle of 6-70° can be achieved, and the 
angle of the measurement arm is independent of the incidence angle. This instrument is able to characterize a sam- 
ple's complete polarization properties in the form of a Mueller matrix by impinging light on it with any state of (fully 
polarized) elliptical polarization. The Stokes vector of the emerging light is measured by a complete polarization state 
analyzer. The sample's Mueller matrix is obtained after complete polarization measurements for at least four different 
input polarization states. The Stokes coordinate system is chosen such that +Q corresponds to vertical and horizontal 
linear polarization, respectively. Positive Stokes U is oriented at 45° clockwise from +Q when looking along the 
beam. -t-V is oriented counterclockwise. The polarization state generator (PSG) consists of a vertical polarizer (0°), 
followed by a compensated liquid crystal variable retarder (LCVR) from Meadowlark Optics at 27.4° (clockwise as 
seen looking downstream) and one at 72.4°. Two LCVRs, the first at -72.4° and the second at -27.4°, followed by a 
vertical polarizer, form the polarization state analyzer (PSA), which is indeed a mirror image of the PSG. For every 
individual LCVR, the retardance curve as a function of voltage is unambiguously determined by measuring the trans- 
mission with the LCVR between parallel and crossed polarizers, with the fast axis at 45° with respect to the polarizers, 
sweeping the voltage with steps of 50 mV. The modulation / demodulation scheme is given by the 16 retardance 
combinations of 1 /8/1 a nd 5/8/i at the LCVRs. T his set-up minimizes the error propagation to the determined Mueller 



matrix, as explained bv lDe Martino et al.l (120031) . To reduce the i nfluence of systematic polarimetric errors due to im- 



perfect components, the eigenvalue calibration method (ECM) bv lCompain et al. I (ll999l) is implemented, which is the 



most accurate method currently known in the literature. The ECM implies the determination of the system matrices, 
describing the PSG and the PSA, on the basis of measurements of a few reference samples with Mueller matrices 
with known parameter dependencies. The reference samples are a Glan Thompson polarizer at 28° and 73°, and a 
Soleil Babinet compensator with A/4 retardance at 23° and 68°, both used in transmission. This choice is based on 
simulations of the calibration, showing that mutual angles of 45° yield the most accurate results, and these exact an- 
gles maximize the accuracy by producing maximum output intensities, making the set-up least sensitive to noise and 
detector non-linearity. In front of the PSG, a broadband halogen lamp, a coUimating achromatic lens and a diaphragm 
produce a collimated beam, 1 cm in diameter, centered on the sample. At the end of the PSA, the light is focused onto 
a photodiode detector. A similar detector directly measures the lamp intensity at the PSG, such that lamp intensity 
fluctuations can be eliminated by normalizing all measurements of the PSA's detector to the reference detector. Differ- 
ent interference filters are used to investigate the wavelength dependence of a sample's Mueller matrix. The relatively 
low photon flux of the light source limited the choice of interference filters for this research to /l=500, 550, 600 and 
650 nm, all with 10 nm band-pass. Since we are primarily interested in the polarization properties of a sample and 
not in its photometric properties, all measured Mueller matrices are normalized to their [1,1] elements that only give a 
measure for its transmission or reflectivity. With this set-up, a polarimetric sensitivity (i.e. the random noise level) of 
2 ■ 10 is achieved in the normalized Mueller matrix elements after a measurement taking 2 seconds per modulation 
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/ demodulation setting, at one incidence angle and wavelength. The accuracy of the determined Mueller matrix per 
measurement is determined by the root mean square (RMS) of the elements that have a perfectly known theoretical 
value of (and if possible also by the RMS of the deviations of the elements that are theoretically identical to the [1,1] 
element, i.e. normalized 1). All measurements presented here exhibit an absolute accuracy of ~ 1% of element [1,1], 
i.e. ±0.01 per normalized Mueller matrix element, and the errors are mainly systematic as is shown in Fig.[T]3. This 
is likely due to the temperature sensitivity of the insulated LCVRs that are measured to be stable to within 0.25°C 
during measurements and calibrations, which corresponds to a maximum retardance variation of 0.625 nm according 
to the supplier's specs. As a consequence, the elements of the normalized system matrices vary up to 0.01. The error 
propagation to the measured Mueller matrix fully depends on the Mueller matrix of the sample and the calibration 
components. Note that fixed retardance deviations, due to limited LCVR calibration accuracy, are detected by the 
ECM. 

220 ±10 nm of aluminum were evaporated onto a 5 x 5 cm square glass substrate under vacuum conditions 
(1.33 ■ 10"^ Pa). The aluminum thickness is determined from step height measurements, performed at different lo- 
cations on the mirror Hence it follows that the uncertainty in the thickness is a direct measure of the flatness of the 
surface. The beam is assumed to be large enough to eliminate the effect of local roughness. The size of the footprint 
changes as a function of incidence angle, but the mirror's macroscopic properties are also assumed to be uniform. The 
characterization of the aluminum oxide layer was consecutively performed by means of reflection ellipsometry. At five 
different intervals after evaporation, and at four wavelengths (see Table [T]for observation and calibration details), the 
Mueller matrices of reflection off the aluminum mirror at 14 angles of incidence (see Fig.[Tll were measured. Both be- 
fore and after the measurements at one wavelength, the ECM was performed to determine the system matrices. Fig.[Tt 
shows the normalized Mueller matrices at 768.8 hours after evaporation, at 600 nm, with the calibration performed 
after the measurements. The top right and bottom left squares, together with element [2,2], determine the accuracy of 
the measurements represented as lo" error bars for all measured Mueller matrix elements, as described above. Table[T] 
shows all accuracies and the ellipsometer's stability over time, both long and short term, and over wavelengths. 
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Fig. 1. — a. Mueller matrices of reflection off an aluminum mirror 768.8 hours after evaporation, at 600 nm. The solid 
curves are the least squares fits to a Mueller matrix model of reflection ofl" aluminum with (black) and without (grey) 
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Table 1: Absolute accuracy of the measured normalized Mueller matrices in percentages of element [1,1]. At every 
wavelength, the values with the calibration before and after the measurements are shown, respectively. 





Absolute accuracy {%) 


Time 




Wavelength [nm +5] 






[hours ±0.3] 


500 


550 


600 


650 


1.3 


1.16 0.88 


1.38 1.33 


1.15 1.40 


1.15 


1.13 


22.0 


1.39 0.95 


1.50 1.18 


1.36 1.18 


1.47 


1.42 


69.5 


0.87 0.82 


1.13 1.02 


1.07 0.92 


1.00 


0.80 


244.2 


1.03 0.96 


1.28 1.09 


1.11 1.01 


1.02 


1.03 


768.7 


0.80 0.77 


0.90 0.79 


0.88 0.75 


0.92 


0.79 
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3. Modeling 

The measurements are interpreted with a theoretical Mueller matrix of bulk aluminum with a thin aluminum oxide 
layer on top of it. This model consists of t he usual reflection Mueller matrix an d complex reflection coefficients (see 
Kelleiil2002l) from thin-film matrix theory (" MacleodllTiei iBorn & Woljll975 l). It considers a bulk metal (subscript 
b) with an index of refraction ni,(A) = h{A) - ik(A) with k{A) > 0. This minus sign follows from the fact that in the 
above thin-film calculations the electric field is proportional to exp[ia>(t - ^z)] instead of exp[/w(-f + 'jz)]- Hence 
it follows that a wrong sign makes k a growing factor, while it should be a damping factor, because of conservation 
of energy. On top of it a dielectric, amorphous thin film (subscript /) with thickness d f is located. The surrounding 
medium is represented by subscript m, which in our case is air with «„, = 1 .00. The values of h(A) and the refractive 
indices of the aluminum oxide fil m (rif) at the applied wavel engths are shown in Table|2l The values of k(A) are widely 
varying throughout the literature ( Harrington & Kuhnl2008 ) and are therefore assumed unknown and fitted along with 
the AI2O3 thickness df. 

The reflection coefficients of p- and s-polarized light, as a function of the angle of incidence Oo, are given by 



(1) 



where 



E,„ - cos (5/ + /2j sin (5/, 
H„, - rji, cos 5/ + irjf sin 5/, 



(2) 



with 



2n 

6f — — nf(A)df cos Of. 
A 



(3) 



The angles 0/ and 6b can be calculated from the angle of incidence using Snell's law. 



sin 6*0 = n f(A) sin 0/, n&(/l) sin 0b - n f{A) sin 0f. 



(4) 



In the above formulae. 



Tlbjjn - nbj,m{^) COS 0bJ,O, 



(5) 



for p polarization, and 



J]b,f,m 



nb,f,mi'l) 



COS 0b,ffl 

for s polarization, respectively. The resultant normalized reflection Mueller matrix is given by 



(6) 
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Mr(«„,(^), nb(A),nj(A),df, 0,,, A) = 



2yR„R,cos(e„-e,) 



1 







R,,+«,, 

1 












R„+R, 
2 '\/RpR^ sin(e/,— e,) 






2 -URpR, cos(e/,-e,) 



(7) 



where 



=arg(rp,,). 



Weighted least squares fits are performed to determine the aluminum oxide thickness and the wavelength-dependent 
values of the imaginary part of the refractive index of aluminum for each time period. Because n j is assumed to be 
sufficiently well known, and an extra degree of freedom causes unstable and unrealistic fit values, «/ is not included in 
the fit (see Table|2]for the literature values). The sum of the quadratic differences between every determined Mueller 
matrix element and the corresponding element in the model is minimized with respect to df, kb{A - 500), kb{A - 550), 
kbiA - 600) and kb{A - 650). The inverse of the accuracy per calibration determines the weight of the corresponding 
Mueller matrices in the fit. As the errors of the determined values are difficult to establish analytically, Monte Carlo 
simulations are performed to determine the spread in the obtained d and k. Each Mueller matrix element (except ele- 
ment [1,1]) at each wavelength and incidence angle is randomly either increased or decreased with the accuracy value, 
after which the fit is performed again. This is repeated 1000 times, and the standard deviations of the determined 
parameter values are used as Icr accuracies of d and k, even though the measurement errors are not Gaussian. 

It has been shown th at the presence of a thin dielectric layer causes conside rable changes in the polarization 



properties of a mirror (e.g. lBurge & Bennetllll964tlSankarasubramanian et alJll999 ). If the layer is neglected, the fit 



results in angle dependent pseudo-indices of refraction n^(A) = hP{A) - ik^iA) that partly correct for the incomplete 
model. To quantify these indices of this mirror, the fit and Monte Carlo simulations are also performed with the 
thickness set to 0, with hP(A) and A:''(/l) as free parameters. Furthermore, the fit is performed with the imaginary part of 
the refractive index of aluminum as the only fit parameter (denoted by kPP{A)), in case one adopts the literature values 
of n, and the aluminum oxide layer is neglected. 

Summarizing, the three types of performed weighted least squares fits are: 
Fit 1 . fit to (i and k, literature values of h 

Fit 2. fit to h and k (denoted by hP and k'', respectively), no aluminum oxide layer in model 
Fit 3. fit to A: (denoted by k''''), literature values of «, no aluminum oxide layer in model. 
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Fig. 2. — Result of Monte Carlo simulations of a 768.8 hours old mirror: 1000 fits of the data with random variations 
in d and k with an amplitude equal to the measurement accuracy (see Table [TJ. The standard deviations of d and k 
determine their accuracies. 



-9- 



4. Results and Discussion 

In Fig. [T^ one set of Mueller matrices is shown as well as the best fit for d and k for that day, with and without 
aluminum oxide layer in the model. The residuals (Fig.lTJ)) show that with the literature value of the real part of the 
index of refraction of aluminum, with a pseudo-value of the imaginary part of the refractive index, the extended model 
is required to accurately describe elements [1,2] and [2,1]. The result of Monte Carlo simulations of the measurements 
of one day is shown in Fig. ID Despite the fact that the measurement errors are not Gaussian, the obtained distributions 
of d and k are inspected to have a distinct Gaussian shape. It shows the error propagation to d and k, and it shows that 
the parameters are correlated, but not to the level that a revision of k can compensate for the aluminum oxide layer: 
its thickness d is significantly larger than zero. Fig. [3]shows the determined values of k(A) and kPP(X), as well as their 
weighted averages per wavelength. It is clear that for Fit 3 (not taking the aluminum oxide layer into account) the 
determined values for k significantly decrease with time, whereas the complete model (Fit 1) yields constant values 
for k. Table [2l shows the weigh ted averages of A:(/l) and the pseudo-index values, and h''(A). In agreement with 



Sankarasubramanian et al.1 (Il999h and lJoos et al.1 (120081) . the absolute values of all pseudo-indices are smaller than the 



actual values. 

The determined aluminum oxide layer thicknesses at different times after evaporation are shown in Fig. |4] with 
linear and logarithmic axes for the time after evaporation. The thicknesses are the best fit values for each measurement 
day, and the standard deviations therein are determined as the standard deviations of the simulated thicknesses per 
day, i.e. the standard deviation of d at one wavelength in Fig. |2] for each measurement day. As the already present 



alumi num oxide prevents oxygen from reaching the aluminum, the layer growth is likely a logarithmic process (iMott 



19391) starting out from small AI2O3 seeds. In order to describe the AI2O3 layer thickness as a function of time, a least 
squares fit is made to the function d(t) - a + blnit). The result, a - 3.74 + 0.18 nm and b - (6.6 ± 3.6) ■ 10"^, is 
overplotted in Fig. H) The corresponding reduced chi squared is - 1.26. As the growth has obviously reached a 
Umiting layer thickness after the first day, the weighted average is also calculated without the first measurement point. 
The result for the long term layer thickness is c/ = 4.12 + 0.08 nm, with^^ = 0.56. 

The ten times larger thicknesses found by Sankarasubramanian et al. I (Il999l) can possibly be explained by the 



applied values of the refractive indices of aluminum and aluminum oxide, including the sign of the imaginary part of 
the index of aluminum. In fact, by applying the opposite sign of k{X), we obtained thicknesses of about 50 nm when 
fitting our data. 
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Table 2: Left: literature values of refractive index of amorphous aluminum oxide («/) and real part of refractive index 
of aluminum (n,) at the applied wavelengths (A). 

Right: real and imaginary parts of refractive index of aluminum, determined with three methods (see Section [3]i: Fit 
1 : a fit of the measurements to k and d (literature values of n), Fit 2: to «p and A;P (no aluminum oxide layer in model), 
and Fit 3: to fcPP (literature values of n, no aluminum oxide layer in model). 





Literature 


Determined 


A [nm] 




k 


«P fcP 


fcPP 


500 
550 
600 
650 


1.61 0.769 
1.61 0.958 
1.60 1.200 
1.60 1.470 


5.88 + 0.02 
6.30 + 0.03 
6.85 + 0.03 
7.33 + 0.03 


0.574 + 0.008 4.93 + 0.01 
0.675 + 0.010 5.32 + 0.01 
0.848 + 0.011 5.80 + 0.01 
1.07 +0.013 6.23 + 0.01 


4.94 + 0.01 
5.33 + 0.01 
5.82 + 0.01 
6.25 + 0.01 



f +0.01, 
* +0.01, 



Eriksson et al.1 (119811) 

Lidd (l2008b . Linearly interpolated 
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Fig. 3. — Fits for k (black circles) and k^i^ (grey squares) at different wavelengths with and without aluminum oxide 
layer in the model, respectively. The lines are the weighted means per wavelength for k and k^P (dashed) (see Table|2|. 



200 400 600 800 

time after evaporation [hours] 

5.0 
4.5 



4.0 
3.5 



1.3 22.0 69.5 244.2 768.7 9000 

time after evaporation [hours] 

Fig. 4. — Results for the AI2O3 layer thickness as a function of time (linear and semi-logarithmic, respectively). The 
solid curve is the least squares fit to a logarithmic function, and the shaded region shows the error bars of the fit. The 
dashed Unes show the weighted average and its error bars of the last four measurement points. 
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5. Conclusions 



We presented the results of Mueller matrix ellipsometry of reflection off an aluminum mirror. Comparison with 
a thin film model shows the necessity of considering the amorphous aluminum oxide layer to accurately describe 
elements [ 1 ,2] and [2, 1 ] of the Mueller matrix. These elements are very important in astronomical polarimetry, because 
they describe the production of linear polarization (i.e. instrumental polarization). Without the aluminum oxide in the 
model, this error is about 0.02 at a wavelength of 600 nm and an incidence angle of 45° (see Fig. [T]|. The cross-talk 
components t/ <-» V (elements [3,4] and [4,3]) can be represented by both the model with and without the oxide layer, 
because both d (aluminum oxide layer thickness) and k (imaginary part of the refractive index of aluminum) can affect 
the mirror's retardance. This is also the reason why d and k are found to slightly correlate in Fig. [21 the presence of the 
layer can to some degree be accounted for by a decrease of k. This is also observed in Fig. [3] This has an interesting 
application in the sense that a mirro r's retardance , and to some extent the induced polarization, can be modified by 
overcoating it with a dielectric layer (lKellei]l2002l) . 



The determined layer thicknesses at 5 different days after evaporation indicate logarithmic growth in the first few 
h ours af t er eva poration, but it remains stable o n the long-term at ~ 4 nm, which is consistent with the theoretical model 
of lMottI (119391) and the experimental results of lJeurgens et alj (l2002h . 



With the obtained values of k(A) (see Table |2|, d - 4.12 + 0.08 nm and the presented formulae (equations ([B- 
([8]l), a normalized Mueller matrix model of a realistic aluminum mirror can be constructed with an absolute accuracy 
of i.e. ±0.01 per normalized Mueller matrix element. Our measurements were performed on a clean mirror, 
disregarding the influence of dust that is always present on telescope mirrors. The polarization properties of a dusty 
aluminum mirror will be presented in a forthcoming paper (Snik et al. 2009; in preparation). 



The authors thank the reviewer, David Harrington, for his constructive and detailed comments on the manuscript. 
We also thank Jatin Rath for evaporating the mirror. 
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